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1—5 

I We present a simple, yet accurate approximation for calculating the cosmic 

microwave background anisotropy power spectrum in adiabatic models. It 
^ , consists of solving for the evolution of a two-fluid model until the epoch 

O ! of recombination and then integrating over the sources to obtain the CMB 

^ I anisotropy power spectrum. The approximation is useful both for a physical 

^ I understanding of CMB anisotropies, as well as for a quantitative analysis 

I of cosmological models. Comparison with exact calculations shows that the 

^ I accuracy is typically better than 20 percent over a large range of angles 

'-^ ' and cosmological models, including those with curvature and cosmological 

^ I constant. Using this approximation we investigate the dependence of the CMB 

?-j I anisotropies on the cosmological parameters. We identify six dimensionless 

^ I parameters that uniquely determine the anisotropy power spectrum within 

our approximation. CMB experiments on different angular scales could in 
^ ^ principle provide information on all these parameters. In particular, mapping of 

' the Doppler peaks would allow an independent determination of baryon mass 

density, matter mass density and Hubble constant. 



Subject headings: cosmology-cosmic microwave background 



1. Introduction 
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Observations of fluctuations in cosmic microwave background (CMB) can provide 
important constraints on cosmological models. Large angular scale (> 10°) observations 
probe the initial conditions, in particular the amplitude and the slope of primordial power 
spectrum ( [Smoot et al. 1992| ; Gorski et al. 19^^ ; [Wright et al. l9I?l ). These scales 
could also provide information on the geometry and the matter content of the universe 
( [Kofman fc Starobinsky 1985| ; [Kamionkowski fc Spergel 1994| ; pugiyama fc Silk 1994|) . 
However, theoretical interpretation of measurements on these scales is complicated by 
cosmic variance and this intrinsically limits the accuracy with which these parameters can 
be estimated using large angular scale measurements alone. While small scale measurements 
suffer less from cosmic variance, their interpretation is complicated by the microphysics 
during recombination and/or reionization. Theoretical models often give wildly different 
predictions for the anisotropy power spectra when the parameter values are only slightly 
changed, while some combinations of parameters seem to provide nearly identical spectra 
(e.g. [Bond et al. 19931 ). The purpose of this Letter is to clarify which combinations of 
physical parameters affect the CMB fluctuations and what are the physical processes that 
lead to these fluctuations. Our two-fluid model for adiabatic fluctuations, presented in 
§2, generalizes previous theoretical approximations that investigated CMB fluctuations in 
the limiting cases of large and small angles (^achs fc Wolfe 1966 ; J0rgensen et al. 1994[) . 
The model is accurate enough that it can be used for a quantitative analysis of various 
models, yet it is also simple enough that it can clearly separate between different physical 
processes that affect the CMB fluctuations. In §3 we use this model to identify the physical 
parameters that can be determined using CMB measurements over a large angular range, 
extending previous studies that were limited to a smaller range of angles and/or parameters 
( [Bond et al. 199"4t [Kamionkowski, Spergel fc Sugiyama 199^1 ; [Gorski fc Stompor 1994[ ; 



^ugiyama fc Silk 1994[; [Muciaccia et al. 1993[). In §4 we present the conclusions. 



2. Method 



In this section we present a method for computing the CMB anisotropies that is a 
generalization of analytic methods flrst introduced by Sachs & Wolfe (1966) and is based 
on the line-of-sight integration along the photon past light cone. We will assume that the 
photons and baryons are tightly coupled prior to recombination, which will allow a simple 
two- fluid description of perturbations ([Peebles fc Yu 1970]) . Our analysis will be restricted 
to the linear perturbation theory of adiabatic perturbations and we will neglect any possible 
tensor contributions. Different theoretical models will be compared using angular power 
spectrum of CMB anisotropies, C; = (| aim P), where aim is the multipole of temperature 
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anisotropy A(n) = J2i,ni(^imyim{n) and Yim{n) is the spherical harmonic. Expected CMB 
anisotropy for a given experiment can be calculated from (A^) = I]«>2(2^ + '^)WiCi/4:H, 
where Wi is the experiment window function. 



The temperature fluctuation A(n) in the direction n can be expressed as a line-of-sight 
integral, 

/■TO A 

A(n)=/ [il{^+^ + n-v-l)+2^]e->'dT. (1) 
Jo 4 

Here r is the conformal time with the value Tq today, is the gravitational potential, 
6^ is the photon density perturbation and is the electron velocity. We introduced the 
Thomson opacity along the past light cone /i(r) = /i(r')(ir' with fi = aXeUeCTT, where 
a is the expansion factor, the ionization fraction, ng the electron number density and 
(Jt the Thomson cross section. The above expression is written in the gauge invariant 
formalism ( [Kodama fc Sasaki 19841 ). We neglected the anisotropic stress contribution and 



the unobservable monopole contribution arising from the local gravitational potential. In 
the limit of infinitely thin LSS the function fie~^ approaches to a Dirac delta-function 
6{t — Tree), whcrc Tree dcuotcs the couformal time at recombination. Equation |l| then 
reduces to (Kodama & Sasaki 1984) 

A(n) = (j){rree) + ^-^^ + u ■ v,{Trec) + 2 T 0(r)cir. (2) 



The velocity term in equation ^ can be rewritten using the linear theory approximation 
Vh = Vipb into n ■ Vb = dipb/dr, where r is the radial coordinate. We may decompose the 
sources in equation ^ into an orthonormal spherical basis set with Fourier amplitudes 0(fc), 
6^{k) and Vb{k), where Vb{k) = iipb{k)/k. After the angular and ensemble averaging we 
obtain the following expression for the multipole moments ( [Kodama fc Sasaki 1984|) , 



oo 

2 TDiv.\rr/u\ n2. 



Q = Att k'P{k)T{k)Dfdk 
Jo 

A = (0 + -j)jl{kTQ - kTree) + vj[{kTQ - kTree) + 2 / dTji{kTQ - kT)F{T), (3) 

where ji is the spherical Bessel functions and ][ its derivative. All the perturbed quantities 
are evaluated in fc-space at T^ee- P{k) denotes the primordial power spectrum of potential 
0, usually expressed as a power law P{k) oc For later purpose we introduced the 

function T{k), which incorporates the damping effects. In the case of a non-fiat universe 
the functions ji need to be substituted by their appropriate generalizations. This is 
unimportant for / > ^|1 — ^^o|^^^ ^^nd the results presented here will be valid for non-flat 
universes, provided that the relation between angles and physical sizes is expressed using 
angular size distances. Last term in equation ^ gives the so-called integrated Sachs- Wolfe 
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(ISW) contribution. The time dependence of the potential is denoted with F{t), where 
-F(r) = 0(r)/0(ro). This term vanishes in a flat, matter dominated = 1 universe, but is 
present in the case of a vacuum energy dominated universe ( [Kofman fc Starobinsky 1985| ), 
curvature dominated universe ( Kamionkowski fc Spergel 1994| ) or when the universe is in 
transition epoch from being radiation to being matter dominated (|Kodama fc Sasaki 1986| ). 

To calculate the anisotropy power spectrum we need to evaluate the source contributions 
in equation ^ at the epoch of recombination. The photon evolution equations in fc-space are 
given by (|Ma fc Bertschinger 1994] ; [Wilson fc Silk 1981| ; [Bond fc Efstathiou 198|) 

4 • k6 

6^ = --kv^ + 4:(p, = + fi{vb - v^) + kcj). (4) 

We also need the evolution equations for baryon and CDM perturbations, 

(jfo = -kvi, + Vh = --Vb + TT-^l^iv^ - Vb) + k(f) 

a Spb 

5c = —kvc + 30, Vc = Vc + k(j), (5) 

a 

where pb and are the baryon and photon mean densities, respectively. 

The energy and momentum constraint equations give the equations for and 0, 

47rGa2 3a/ • d A-nGa'' f 

-(P+^), = — 0+ , (6) 



k"^ ak a k 

where p = (P7 + Pv)5.y + pbh + pc^c and / = ^{p-y + Pv)v^ + pbVb + pcVc- Here p^, and pc 
are the neutrino and CDM mean densities, respectively. We replaced neutrino density and 
velocity perturbations with the corresponding photon perturbations. This becomes invalid 
on small scales due to the free-streaming of neutrinos, but does not affect significantly the 
final results. We also neglected the anisotropic shear and the possible curvature terms. 

The above equations are supplemented by the Friedmann equation, which at early 
times (when a possible cosmological constant or curvature term can be neglected) is given 
by 

,d,2 87rG'a^, , 

{-) =^^{P1 + Pu + Pb + Pc). (7) 

The solution to this equation is 

y = — = (axY + 2ax, x = ( — = — , (8) 



Oj^q \(Xrpr/ 



where aeq = {pj + Pu)/{Pb + Pc) ~ 4.2 x 10 ^^^h ^ (assuming three flavors of massless 
neutrinos), a~^^ ^ 1100 for the standard recombination, = arec/deq, = 0,b + is 
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the value of matter density today in units of critical density and h is the value of Hubble 
constant today in units of lOOkm/s/Mpc. 

We will now assume the tight coupling limit /i ^ 1, which is a good approximation on 
scales larger than the Silk damping scale ( pilk 196"8| , Peebles fc Yu 19701) . In this case the 
photons and baryons are coupled into a single fluid with 6b = ^S-y and Vb = v^. The above 
equations rewritten in terms of dimensionless time x and dimensionelss wavevector n = kTr 
become 



Sr. 



-KVc + 30 . 



2 hhV^ ~t~ 



l+y 



-rjVc + n4> 



+ Vh) 



-vvbVy + — + f^n^ + Vb) 



0= -\{ri/ kY{5 + ?>riv/ k) , = -77^ + 



?)rfv 
2k 



v-yi^ + y - Vc) + VcVc 



l + y 



(9) 



where the derivatives are taken with respect to x, yb = 



3x7 / 
8 ^ 



1.68^», 



yc = ^y = (1 — ^/1.68)?/ and r] = 2a{ax + + 2ax). Equations ^ are a coupled 

system of 4 first order differential equations.^ The appropriate initial conditions at a; << 1 
(when the universe is radiation dominated) and kt] << 1 (when the mode is larger than the 
Hubble sphere radius) are 



6^ 



-20(1 + 



16' 



6c 



V 



5^^2^:2(1 + ^)0 



V(| + y) 



(10) 



The above equations need to be evolved until Xrec = [(«^ + 1)^^^ — l]/a- 
temperature anisotropy expressed with the dimensionless variables is given by 



The 



Ci = 471 A / K"T(K)dlnK 
Jo 



+ ^ + 2A0)j j,(/txo) 



VjJi[KXo} 



3^n— 1 



where xq is the angular distance to the LSS in units of and we assumed P{k) = Ak 
The term A0 = [2 — 8/y{xrec) + lGxrec/y^{xrec)]/10y{xrec) arises from the ISW effect due to 
the potential varying with time during the transition period from the radiation dominated 



^ In actual numerical implementation of these equations we find that for a stable numerical integration 
it is better to compute (f> using its time evolution in equation 0, rather than computing it from the sources. 
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to the matter dominated universe ([Kodama fc Sasaki 1986 ). For simplicity we dropped the 



ISW contribution from possible curvature or cosmological constant, which is only important 
at the lowest values of / (/ < 10). 

The damping transfer function T(/t) is approximately unity for low values of Z (/ < 200), 
but gradually decreases afterwards. Its main contributions come from the Silk damping and 
from the finite width of LSS. The first effect can be calculated analytically by expanding 
equations ^ and |^ to second order in /i and neglecting the effects of gravity and expansion. 
In the matter-dominated era one obtains T(k) oc exj9(— 2k^x^) ( [Fugukita et al. 1990|) , where 
Xs is the Silk damping scale in units of r^, Xs = 0.6fi^/^l]6"^/^a3/4/i-i/2. Second effect can 
be analytically estimated by performing the line-of-sight integral in equation ^ in the limit 
where the sources are slowly changing over the timescale on which the visibility function 
fie~^ is non-negligible ( |j0rgensen et al. 1994| ). Visibility function can be approximated as 
a gaussian {2Tca'^x'^^J~^^'^exp[—{x — XrecY /'^{(^Xrecf], where for standard recombination 
a 0.03.0 In the limit k^qXq » I (where is the wavevector which gives the dominant 
contribution to Q), we obtain the damping factor exp{—H^a'^x1^^. Therefore, the damping 
effects can be written as 

T{k) ^ e-'^'(2^'+"'^-). (12) 

This works reasonably well for the standard ionization history. Note however that if the 
limit k^qXq >> / is not satisfied, then the damping due to the finite thickness of LSS is not 
exponential, but is proportional to /t"^. This will be the case, for example, in reionized 
models. 

Equations ^|T2] are all is needed to evaluate the temperature fluctuations. Although 
equations ^ cannot be solved analytically in general, they have analytic solutions in the 
limits of small and large k. In the first limit where the modes are larger than the Hubble 
sphere radius the amplitude of perturbations at a given time is a constant (figure p. This 
gives the standard Sachs- Wolfe expression for CMB fluctuations, as can be verified by 



evolving the initial conditions in equations 10 into the matter-dominated era and neglecting 



the velocity term in equation |Tl|. In the second limit (large k) the equations can be solved 
using the WKB approximation and the solution is given by the acoustic oscillations of 
photon-baryon plasma ( |J0rgensen et al. 1994| , [Padmanabhan 19931 ). In the intermediate 



regime, which is of main interest for us, the equations need to be solved numerically, but 
the physics can be well understood by the two limits above. As shown in figure ^ equations 
^ give an excellent approximation to the exact results over a large range of wavevector k. 



•^Both a and arec are weakly dependent on cosmological parameters. Moreover, the ISW visibility function 
differs from the Thomson scattering visibility function, which leads to a different damping of ISW term. Both 
effects will be neglected here. 



- 7- 



Although the above system of equations is particularly useful for the standard 
recombination scenarios, one can also use it to calculate anisotropy power spectrum in 
reionized models. As one can see from equation |I|, the primary fluctuations will be 
suppressed by a factor exp[—fi{xrec)], where fi{xrec) is Thomson opacity at recombination. 
In addition, there will be secondary fluctuations generated at the new last-scattering 
surface, which can be calculated using the same method as above, except that one needs to 
replace with Vc in the regime where the Compton drag is negligible. Since the modes 
larger than the Hubble sphere radius at the new LSS do not evolve in the matter-dominated 
regime this simply regenerates the Sachs- Wolfe expression for low values of 1. On smaller 
scales the thick new LSS damps the secondary fluctuations and in many scenarios these 
become negligible (although on arcminute scales the second-order terms may become 
important, [Vishniac 19871 ). Provided that one is interested in degree angular scales, then 
the effect of reionization is to suppress the fluctuation power spectrum relative to the large 
scales by a factor of exp[—2fj,{xrec)]- 



3. Results 



The comparison between our approximation and the exact solutions of perturbed 
Boltzmann equation (Seljak & Bertschinger 1994; Seljak 1994) is presented in figure ^ 
for several cosmological models. All of the multipole moments are normalized relative 
to Cio, which is approximately fixed by the COBE experiment and where the curvature 
effects and ISW effects due to Qm 7^ 1 can be neglected. One can see that the agreement 
is excellent over a large range of I. The deviations at large / arise because of improper 
treatment of damping effects, while the deviations at small I can be attributed to the 
neglection of neutrino anisotropic shear. Another effect that introduces small deviations is 
the dependence of a^ec on cosmological parameters, which slightly offsets the position of 
the peaks. Nevertheless, our approximation correctly predicts the positions and amplitudes 
of Doppler peaks (also called acoustic or Sakharov oscillations) with a 20 % accuracy over 
most of parameter range. 

Given the high accuracy of our model we may now investigate how the anisotropy 
power spectrum depends on the cosmological model. Our goal is to identify the parameters 
that can be determined using the CMB measurements and to give a physical understanding 
of how they affect the anisotropics. In equations ^ the free parameters are a = 
and ^. We can replace them with physically more relevant parameters Vtmh"^ = (a/21. 5)^ 
and fife/i^ = (a/21. 5)^^. In addition to these two we have the parameters /i(xrec), xq 
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and Xg . These six parameters will uniquely determine all the CMB power spectra within our 
approximation. Of the parameters above, n, fi{xrec) and Xs all suppress the power on small 
scales relative to large scales (for n < 1). The suppression is different in the three cases, 
being proportional to k""^ (neglecting the possible tensor contribution), exp[—fi{xrec)] and 
exp(— respectively. This, in principle, allows to separate the different suppression 
effects and to determine the three parameters separately (but see Pond et al. 199^) . 
Note that these parameters do not change the positions of the Doppler peaks, only their 
amplitude. Since the effects of these parameters are physically transparent, we will restrict 
in the following to the case of n = 1 and fi{xrec) = 0. Moreover, Silk damping is important 
only for large values of / and it can be neglected if one concentrates on the first few Doppler 
peaks. We are thus left with xq, Q^h'^ and flbh'^, which uniquely determine positions of the 
Doppler peaks. 

Position of the first Doppler peak is determined by the angular size of the Hubble 
sphere radius at decoupling, which, expressed in terms of our variables is given by 
(1 + Q;~^)^/^xo. Since a >> 1 for typical values of Vlmh?^ position of the first Doppler peak 
mainly depends on xq. Assuming Xrec « Xq we have Xq = (fimC^rec)"^^^ for the model 
with negligible cosmological constant and xq = il'^^a~^J'^ for the model with neghgible 
curvature. The latter result shows that the first Doppler peak only weakly depends on 
fi^ = 1 — ^m- This is because the angular size distance at large redshifts scales with in 
approximately the same way as does the Hubble sphere radius at decoupling ( |Vittorio~S^ 



|Silk 1985| ; |G6rski fc Stompor 1994| ). This is, however, not true in general and the position 
of the first Doppler peak depends on Qm and Qx when both curvature and cosmological 
constant are important. The value of I at which the maximum of the first Doppler peak 
lies is given approximately by 6xq. Figure ^ja compares Q's of curvature and cosmological 
constant dominated models with those of = 1 model at fixed values of flmh"^ and Qfyh"^. 
One can see that the position of the first Doppler peak can accurately determine flrn in 
curvature dominated universe ( [Kamionkowski, Spergel fc Sugiyama 1994| ), but cannot 



precisely determine Qm in cosmological constant dominated universe ( [Bond et al. 199^ ; 
|G6rski fc Stompor 1994| ). However, even in this model the positions of secondary Doppler 
peaks are already significantly displaced relative to each other when Q^n changes from 0.25 
to 1. This would thus allow independent determination of Qm even in cosmological constant 
dominated universe, once Qmh"^ and Qbh"^ are known (see below). 

The dependence of the Doppler peak positions and amplitudes on flhh'^ and Qmh'^ 
is more complicated, since both parameters appear in the evolution equations ^ and 
change the properties of acoustic oscillations. Moreover, the two parameters enter into the 
equations differently and have different physical effects: Qbh"^ is related to the properties of 
photon-baryon plasma and determines its effective sound velocity at recombination, whereas 
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^mh"^ is related to the time evolution of the expansion factor, since it determines the epoch 
of matter-radiation equality. This means that one cannot expect the anisotropy spectra to 
remain invariant under a certain combination of the two parameters and both Q^h"^ and 
Qmh"^ are required for a complete description of the Doppler peaks. Figures |^ and ^ 
show how the Doppler peaks change when one of the two parameters is changing while the 
other is held fixed. If one concentrates on the first Doppler peak then it is not possible to 
determine the two parameters simultaneously, since both increasing Qhh'^ and decreasing 
Vtmh? lead to an increase in the first Doppler peak. The physical mechanisms that lead to 
this are different: while an increase in increases the amplitude of the first wave in 

and + (5^/4 (fig. |l|), a decrease in Vtmh"^ also leads to an increased ISW contribution. 
Once the secondary peaks are observed as well, then different effects of the two parameters 
become significant and allow one to determine the two parameters simultaneously (figs. ^ 
b,c,d). Figure |^d shows how changing Vtmh? at a fixed value of ^ affects the multipole 
moments. Again, since the epoch of matter-radiation equality is changing with Qmh'^, one 
does not expect the multipole moments to remain unchanged and both our approximation 
and exact calculations confirm this. Therefore, changing at a fixed and changes 
the anisotropy power spectrum, contrary to some recent claims ( Pond et al. 199^ ; |G6rsk: 

Stompor 19941 ). 



4. Discussion 



The approximation for calculating anisotropy power spectrum presented here is 
a generalization of the Sachs- Wolfe approximation, which itself is only valid on scales 
larger than the Hubble sphere radius at recombination. By modelling the cosmological 
perturbations as a two-component fluid plasma we extended this approach to all scales. The 
approximation is useful both for developing the physical understanding of processes that 
affect CMB fluctuations, as well as for a quantitative prediction of multipole moments for 
various cosmological models. The main approximations used in our model are a two-fluid 
approximation, neglection of anisotropic shear, a simplified treatment of Thomson scattering 
effects, neglection of curvature effects and neglection of vector and tensor contributions. 
None of these assumptions is essential for the method and one can generalize the approach 
presented here to obtain exact results (Seljak 1994). This will lead to a computationally 
more demanding system of equations, but the main physical effects that lead to the creation 
of Doppler peaks will still be determined by the equations presented in this Letter. 

By rewriting the equations in their dimensionless form we identified all the dimensionless 
parameters that affect the anisotropy power spectra. Measurements of CMB fluctuations 
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cal! only determine these parameters. For example, neutrinos enter into our equations 
indirectly through the Friedmann equation |^ and through the energy-momentum constraint 
equations |^. The presence of a massive neutrino only weakly changes these equations 
and the resultant multipole moments are almost indistinguishable from the ones with the 
massless neutrino. Therefore, the question of whether neutrino has a mass has little hope 
to be answered using the CMB measurements. 

The most interesting aspect of the CMB power spectra is the peculiar pattern of the 
Doppler peaks, which allows a simultaneous determination of Qbh"^ and VLmh?- This would 
provide an independent test of nuclesynthesis prediction of (e.g. [Walker et al. 199"1| ) 



and would also constrain the parameter space on Vt^ and h. In addition, the position 
of the first Doppler peak determines VLm in curvature dominated model. In cosmological 
constant dominated model the position of the first Doppler peak does not allow one to 
determine VL^ accurately, but positions of secondary Doppler peaks could be used to 
constrain VLrn (although for accurate determination exact calculations should be used in 
this case). Another way to break the degeneracy between Vtm-, and h is to determine 
the Silk damping scale x^, which depends only on these three parameters and cannot be 
expressed as a combination of Vti^h? and flmh^- This would require a separation of Silk 
damping from the damping due to the finite thickness of LSS in reionized models and from 
the n < 1 suppression of small scales relative to large scales (including the possible tensor 
contribution). This is possible, because the three effects suppress the small scale power 
differently. Thus, a combination of CMB measurements over a large range of angles could 
be used to separately determine the baryon mass density, matter mass density and the 
Hubble constant value. 

I would like to thank Ed Bertschinger for useful discussions. This work was supported 
by grants NSF AST90-01762 and NASA NAGW-2807. 
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Fig. 1. — Comparison between our approximation (thick lines) and exact solution (thin lines) 
for (solid lines) and + 5^/ A (dashed lines) as a function of k. Silk damping has been 
included according to the expression in the text. Parameter values are 0^ = 0.05, h — 0.5 
and ilm = 1- 
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Fig. 2. — Comparison of anisotropy power spectra between our approximation (thick lines) 
and exact solution (thin lines) as a function of multipole moment I for several different 
cosmological models. The spectra in all figures are normalized relative to Ciq. 
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Fig. 3. — Anisotropy power spectra as a function of multipole moment / for different 
cosmological models. In (a) curvature and cosmological constant dominated models with 
^^?Ti — 0.25 are compared to = 1 model. In (b) Q^h.'^ is fixed at 0.5 and Qbh"^ is varying, 
whereas in (c) flhh'^ is fixed at the nucleosynthesis value and flmh^ is varying. In (d) fli^/fl^ 
is fixed and ^mh^ is varying. In (b), (c) and (d) ^2^^ = 1. In all cases varying the parameter 
changes the pattern of Doppler peaks. 
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